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Abstract 

The Canham-Helfrich free-energy density for a lipid bilayer has drawn considerable attention. Aside from 
the mean and Gaussian curvatures, this free-energy density involves a spontaneous mean-curvature that 
encompasses information regarding the preferred, natural shape of the lipid bilayer. We use a straight- 
forward microphysical argument to derive the Canham-Helfrich free-energy density. Our derivation (i) 
provides a justification for the common assertion that spontaneous curvature originates primarily from 
asymmetry between the leaflets comprising a bilayer and (ii) furnishes expressions for the splay and 
saddle-splay moduli in terms of derivatives of the underlying potential. 
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1 Introduction 

Biomembranes are ubiquitous in nature as they form cell walls. An essential element of a biomembrane 
is a lipid bilayer, which is composed of lipid molecules. These molecules have a hydrophilic head and a 
hydrophobic tail. Due to these properties, when a large number of lipid molecules are placed in a solution 
they self assemble, under suitable conditions, into two-dimensional structures consisting of two leaflets (or 
monolayers) . The lipid molecules are oriented so that the tails of the molecules in each leaflet are in contact 
with each other, while the heads are in contact with the ambient solution; see, for example, Lasic (1988). 
These two-dimensional structures often close to form vesicles and are usually between 50 nanometers and 
tens of micrometers in diameter but only a few nanometers thick, as observed by Luisi & Walade (2000). 
Due to these dimensions, lipid bilayers are usually modeled as surfaces. 

There is a long history of work dealing with specifying a free-energy density for a surface representing a 
lipid bilayer. Over forty years ago, Canham (1970), in an effort to explain the biconcave shape of red blood 
cells, proposed a bending-energy density dependent on the square of the mean curvature. A few years later, 
Hclfrich (1973) also considered an elastic-energy density depending on the curvature of the surface associated 
with a lipid bilayer. In the Canham-Helfrich theory, the free-energy density ip of the lipid bilayer is of the 
form 

i>= X -n(H - H ) 2 + RK, (1.1) 

where H and K are the mean and Gaussian curvatures, n and k are the bending moduli, and H Q is the 
spontaneous mean-curvature. The curvatures H and K are scalar invariants that describe the shape of 



the lipid bilayer. (See Appendix A. 2 for the definitions of H and K.) The splay modulus k, which is 
always positive, characterizes the resistance to changes in mean curvature, while the saddle-splay modulus 
K is related to the resistance to changes in Gaussian curvature. The spontaneous mean-curvature H a is 
determined by the preferred, natural, local shape of the lipid bilayer. Although Canham (1970) and Hclfrich 



(1973) were the first to use a free-energy density of the form (1.1 1 for biomembranes, such expressions were 
considered by Poisson (1812), in the H = case, Germain (1821), and the Cosserat brothers (1909) in 
the study of elastic surfaces. Although not present in the Canham-Helfrich theory, it is also possible to 
introduce a spontaneous Gaussian-curvature K Q , as do Maleki et al. (2012), in which case the last term on 



the right-hand side of (1.1 1 is replaced by k(K — K ). If R and K are both constant, the spontaneous 
Gaussian-curvature can be neglected since the free-energy density is only determined up to an additive 
constant. 



As mentioned above, the spontaneous curvatures H a and K a are determined by the preferred, natural, 
local shape of the lipid bilayer. As Seifert (1997) notes, one of the main causes for spontaneous curvature is 
asymmetry between the two leaflets making up the bilayer. The asymmetry can manifest itself in different 
ways. For example, Dobereiner et al (1999) observe that spontaneous curvature may arise from differences 
in the chemical properties of the aqueous solution on the two sides of the lipid bilayer, while McMahon & 
Gallop (1998) explain that lipid molecules with different configurations may also lead to nonzero H Q and K Q . 

Several derivations of the Canham-Helfrich free-energy density have appeared in the literature. For 
example, in a work concerned with lipid monolayers, Safran (1994) derived ( |1.1| by modeling the lipid 
molecules in a monolayer as springs. In Safran's (1994) derivation, the springs are allowed to change length 
with the bending of the bilayer but may not interact with one another. Granted that the lipid monolayer 
is incompressible, Safran (1994) expressed the length of a generic molecule in terms of the underlying mean 
and Gaussian curvatures. Upon making the additional assumption that the product of the spontaneous 
curvature with the natural length of the molecules is small, he found that the energy governing each spring 
takes a form consistent with ( 1.1 1. 

Paunov et al. (2000) also used statistical methods to derive an energy for a lipid monolayer. Their 
expression for the energy consists of two terms: one associated with bending, which is taken from the work 
of Safran (1994), and the other associated with mixing. The mixing term accounts for interactions between 
the monolayer and the ambient solution. For the particular cases of spherical and cylindrical monolayers, 
Paunov et al. (2000) found that, when expanded up to include only terms up to the square of the inverse 
radius of curvature, their expression for the energy of a monolayer agrees with 

Following Ljunggren & Eriksson (1995), who took the view that the bending energy of a bilayer is due to 
electrical forces, using the Debye-Huckel approximation, Winterhalter & Helfrich (1988) obtained expressions 
for the bending energy of planar, cylindrical, and spherical lipid bilayers. As with the derivation of Paunov 
et al. (2000), Winterhalter & Helfrich (1988) considered only terms up to the square of the inverse radius of 
curvature and, in so doing, obtained expressions for the bending moduli and spontaneous curvature. 

A derivation for the bending energy of a bilayer was also provided by Seifert (1997), who models the 
bilayer as a surface with molecules populating both sides, the molecular density being different on the 
sides. Disregarding the structure of these molecules, Seifert (1997) assumes that the energy of the bilayer 
is determined by the extent to which these densities deviate from an equilibrium density. The densities 
are defined on the midsurface of each monolayer. To obtain an energy density defined on the surface, the 
two molecular densities are projected onto the surface. Since these projections depend on the curvature of 
the surface, the energy density of the surface inherits dependance on the curvature. By assuming that the 
curvature of the surface is small, Seifert (1997) finds an energy density for the surface of the Canham-Helfrich 
variety with no spontaneous curvature. 

In the present paper, we follow an approach used by Keller & Merchant (1991) to derive the energy of 
a liquid surface. To obtain their result, Keller & Merchant (1991) expanded the energy of a mass of liquid 
in powers of the range d of the intermolecular potential. Beyond the first term of this expansion, which is 
volumetric, the second and third terms of this expansion are areal. Of the areal contributions to the energy, 
the first is due to surface tension and is proportional to the area of the underlying surface, and the second 
is proportional to the integral over the surface of H 2 — and, thus, to the bending-energy density arising 
in the plate theory of Foppl (1907) and von Karman (1910) with Poisson's ratio equal to 1/3. 

In our derivation of the free-energy density of a lipid bilayer, we neglect interactions between the lipid 
molecules and the ambient solution and instead focus on interactions between the lipid molecules themselves. 
These molecules are modeled as rigid rods, as in the work of Seguin & Fried (2012). The rods, representing 
the molecules, are distributed over the surface that represents the lipid bilayer. By assuming that the lipid 
molecules are oriented perpendicular to the surface and interact with each other only when they are within 
a distance of d 7 we obtain an expression for the free-energy density for the bilayer. Following the approach 
of Keller & Merchant (1991), we then expand the free-energy density in powers of e :— djl. where £ is a 
characteristic length associated with the bilayer, and identify the Canham Helfrich free-energy density as 
the coefficient of the term of order e 2 . This expansion rests on the assumption that the smallest possible 
curvature of the membrane is much larger than its thickness and, thus, our derivation does not hold for 
membranes that are undergoing crumpling, pitting, or budding. Our derivation also neglects the effect 
of any molecular fluctuations. While these fluctuations are important in crumpling, budding, and pitting 
events, as described by Lipowsky (1990) and Sackmann (1994), all of these events involve the membrane 
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having large curvature and, thus, are not included in our theory. Moreover, we view the Canham-Helfrich 
free-energy density as a mean-field description of a lipid bilayer. In classical mean-field theory, fluctuations 
are neglected. 

This argument provides a direct and easy derivation of the Canham-Helfrich free-energy density. It also 
affords several observations. For one, if the lipid bilayer is symmetric, which occurs when its leaflets are of 
identical composition, and the lipid molecules are uniformly distributed, then the spontaneous curvatures 
H and K Q must vanish. Moreover, since our derivation gives concrete expressions for the bending moduli, 
those expressions might be used to gain insight regarding the sign of the saddle-splay modulus R. 

Our paper is outlined as follows. In Section [2] we describe our approach to modeling lipid molecules, 
write down a free-energy density for a lipid bilayer, and state the result of expanding this energy in powers 
of d/l. The details of this expansion can be found in the Appendices. This is followed, in Section [3] with 
a discussion of the resulting expansion and some consideration of special cases. In Section [4] we discuss 
possible extensions of the method used in Section [2] including some possible approaches to weakening the 
assumptions made there. Finally, in Section [5] we summarize our primary results and make some concluding 
remarks. 

2 Free-energy density of an isolated lipid bilayer 

In this section, guided by the work of Keller & Merchant (1991), we derive an expression for the free- 
energy density of a single closed (lipid) bilayer suspended in solution. Our result is based on the following 
assumptions: 

(i) The thickness of the bilayer is small relative to its average diameter. 

(ii) The (lipid) molecules can be modeled as one-dimensional rigid rods. 

(iii) The molecules do not tilt relative to the orientation of the bilayer. 

(iv) Interactions between the bilayer and the solution are negligible. 

On the basis of Assumption (i), which is commonly imposed in models for closed bilayers (Luisi & Walade 
(2000)), we identify both leaflets of the bilayer with a single orientable surface S. This surface can adopt a 
large variety of shapes; however, being made up of molecules of a finite size, it cannot support arbitrarily 
large curvatures. Let i denote the smallest stable radius of curvature that S is capable of exhibiting. For 
each leaflet i = 1, 2 of the bilayer, we introduce a molecular number density Wi defined on S and measured 
per unit area of S. The total number of molecules in leaflet i = 1, 2 is then given by the integral 

/ Wi(y)da y> (2.1) 
Js 

where da y denotes the area element on S. Taking Wi to be defined on S amounts to assuming that the 
centers of the lipid molecules of both leaflets lie on S, which is consistent with assuming that the bilayer is 
thin relative to its average diameter. In general, the number densities of the leaflets may differ. 

On the basis of Assumption (ii) , the configuration of each such molecule may be described by a point on 
S and a unit- vector- valued director (Seguin & Fried (2012)), with the point representing the center of the 
rod and the director representing the orientation of the rod. Without loss of generality, we suppose that the 
director tips point toward the headgroups of the molecules. We also assume that the interaction between a 
pair of molecules is governed by a potential that depends on the vector between the points and the directors 
and is restricted such that only molecules separated by distances less than some fixed cut-off distance d may 
interact. We assume that the cut-off distance d is small relative to the smallest radius of curvature t the 
bilayer can support, so that d <C £ or, equivalently, 

£■■= j « I- (2-2) 

Given molecules located at two points x and y on S and corresponding directors e(x) and f(y), we 
consider interaction potentials, with dimensions of energy, of the generic form 

<f> d (x - y, e(x),f(y)) = ^{x - y\\ (x - y) ■ e(x), (x - y) ■ f(y),e(x) ■ f (y)), (2.3) 
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Figure 1: Graph of s i-> 0(s, a, 6, c) for fixed a, 6, and c. 



where, to ensure that rods separated by a distance greater than d may not interact, <p must satisfy 



, a, 6, c) = if s > £ for all (a, b, c) 6 



(2.4) 



Notice that $ d depends on d, while does not. That $ d must depend on the vectors x — y, e(x), f(y) 
only through the scalars \x — y\, (x — y) ■ e(x), (x — y) ■ f(x), and e(x) • f(x) is shown by Seguin & Fried 
(2012). This is also a consequence of the requirement that the potential be frame-indifferent. Notice that the 
interaction energy 4>(s, a, b, c) between two molecules can change if the head and tail of one of the molecules 
flips. Interaction potentials of this form may therefore account for differences between the polarities of the 



heads and tails of the molecules. In taking $ d to depend on the cutoff distance d as shown in ( 2.3 ) , we follow 
Keller & Merchant (1991). In addition to a potential $^ that accounts for interactions between molecules 
in each leaflet i — 1, 2, we introduce a potential $f 2 — ^21 that accounts for interactions between molecules 
belonging to different leaflets. Although the specific forms of these potentials may differ, they share the 
same general properties to the extent that they satisfy (2.3 1 and (2.4). A graph of a particular <f> for fixed 



a, b, and c is given in Figure [T] This function is based on the Lennard- Jones 6-12 potential. More tangible 
examples of interaction potentials are provided by Yuan (2010). 

Without loss of generality, we choose to orient S with the unit-normal field that points into the fluid 
adjacent to the head groups of leaflet 1 and denote that field by n. On the basis of Assumption (iii), it then 
follows that the directors of molecules in leaflets 1 and 2 coincide with n and — n, respectively. Bearing in 



mind the cut-off property (2.4), define Sd(x) by 



The integrals 



and 



^11 0*0 = 



S d (x) :={yeS:\x-y\<d}. 



$> d n (x - y, n(x),n(y))W 1 (x)W 1 (y) da, 



S d (x) 



S d (x) 



22 



(x - y, -n(x),-TLi(y))W 2 (x)W 2 (y)day 



(2.5) 



(2.6) 



(2.7) 



then represent the contributions to the value of free-energy density of the bilayer at x due, respectively, to 
interactions between the molecules of leaflets 1 and 2 at a; with all other molecules in their own leaflets. Notice 



that the interaction energy appearing in the integrand of (2.6) gives the interaction energy between a 
lipid molecule at x and a lipid molecule at y. The two number densities are needed to account for the number 
of molecules at x and the number of molecules at y. An analogous explanation explains the presence of the 
two densities in (2.7 1. Similarly, the integrals 



^12(2;) 



S d (x) 



$? 2 (x - y, 11(1), -n(y))W 1 (x)W 2 (y) da v , 



(2. 
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and 



"021 (x) 



$Ux - y, -n(x),n(y))W 2 (x)W 1 (y) da y 



(2.9) 



represent the contributions to the free-energy density of the bilayer at x due, respectively, to interactions 
between the molecules of leaflet 1 at x with all the molecules of leaflet 2 and to interactions between the 
molecules of leaflet 2 at x with all the molecules of leaflet 1. Notice that interactions between the molecules 
in leaflet 1 at x with those in leaflet 2 at x are encompassed by ipi2 (x); analogously, interactions between 
the molecules in leaflet 2 at a; with those in leaflet 1 at x are encompassed by ip2i(x). 

Finally, on the basis of Assumption (iv) , the free-energy density of the bilayer at x is determined com- 
pletely by the sum of the four contributions defined in (2.6 l-(|2.9| and thus has the form 



(2.10) 



On substituting (2.6 |- pi)| into the right-hand side of (2.10), the goal is to expand tp{x) in powers of e 
up to order e 2 , following Keller & Merchant (1991), with the objective of capturing the dependence on the 
curvature. This expansion, the details of which appear in the Appendices, yields 



1 



tp = Vo + ^{H - H o y + R(K - K ), 



(2.11) 



where k, k, H , and K are defined in (A. 50 )— ( A. 53 ), respectively, k and R are of order e 2 and terms of order 
o(e 2 ) are neglected. The quantities <po,K, R, H a , and K depend on the point x in S. Thus, ip mav depends 
on x through not only the mean and Gaussian curvatures but also through the splay and saddle-splay moduli 
and the spontaneous mean and Gaussian curvatures. 



3 Interpretation of results 



The result ( |2.11[ ) deserves some interpretation. The first term on the right-hand side of (2.11 1 only depends 
on the interaction potentials and number densities at x. This term remains in the limit as d (or, equivalently, 
e) approaches zero. The remainder of (2.111 is the Canham-Helfrich free-energy density, augmented by a 
spontaneous Gaussian-curvature term. Hence, (2.111 provides a microphysical derivation of the Canham- 
Helfrich free-energy density for a lipid bilayer. Moreover, the abbreviations ( A. 45 1— ( A. 53 ) describe 
how the spontaneous curvatures H and K and the bending moduli k and R depend on the interaction 
energies between the lipid molecules and the densities of the molecules in the two leaflets. The term ipo 
is independent of the shape of the membrane and not part of the Canham-Helfrich free-energy density. 
However, due to implicit dependence of the number densities on temperature, concentration, and relevant 
electromagnetic fields, that term encompasses effects associated with ambient temperature, concentration, 
and relevant electromagnetic conditions. 

Since the number densities W\ and W2 need not be constant, the spontaneous curvatures H and K 
and the bending moduli n and R generally vary with position x on the bilayer. Since the free-energy density 
can be determined only up to an arbitrary constant, if it happens that both R and K Q are constant, then 
the the term —RK can be dropped without loss of generality. Otherwise, the presence of the spontaneous 
Gaussian-curvature K cannot be neglected. 

To obtain some additional insight regarding the spontaneous curvatures and the bending moduli, suppose 
that the two leaflets are identical, so that there is no leaflet asymmetry. This means that tpu = 1P22 = "012 
and Wi = W%. From (A. 47) and (A. 52 1, we see that in this case C = H D = and, hence, that there is 



no spontaneous mean-curvature. Moreover, if the lipid molecules are evenly distributed on each leaflet so 
that the number densities W\ and W2 are constant, then D = and, hence, by (A. 53), the spontaneous 
Gaussian-curvature K Q also vanishes. These results are in agreement with what other researchers have 
proposed regarding spontaneous curvature. Specifically, Seifert (1997) proposes that spontaneous curvature 
is a measure of the asymmetry of the two leaflets. McMahon & Gallop (1998) suggest that spontaneous 
curvature is a result of the leaflets having different molecular compositions. In the framework presented 
here, possible differences in the compositions of the leaflets are encompassed by allowing tpu and ip 2 2 to 
differ. 
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Ideally, it should be possible to obtain an interaction energy for a given type of lipid molecule and use 



(A. 45 1— ( A. 53 1 to explicitly calculate the spontaneous curvatures and bending moduli for a bilayer consisting 
of these lipid molecules. For example, consider an interaction potential of the form 

1 



(f>ij(s,a,b,c) := -p l3 (s)(c - l)if s > i for all (s, a, b, c) e K X K x R x [— 1, 1], 



(3.1) 



which is taken from the Maier-Saupe (1958) theory of liquid-crystals. Here, pij, i,j = 1,2, are functions of 
the square of the normalized distance s — e~ 2 \x — y\ 2 between the lipid molecules. In this case, it is easily 
confirmed that in^i = 4>ij,3 = and, hence, by (A. 45 1 and (A. 52 1, H Q — 0. Thus, when applied to lipid 



Uij,3 



molecules, the Mai'er-Saupe (1958) interaction potential leads to zero spontaneous mean-curvature. Without 
some concrete information about the functions p^ it is not possible to say anything about the characteristic 
features of the moduli k and R and, thus, in particular, about the sign of R. 
Consider next the simple generalization 



(j}. l3 (s,a,b,c) := ~Pij(s)(c 2 



1 + (a - ay) a + (b- bij) 2 ) 



its>£ for all (s, a, b, c) e R x K x K x [-1, 1] (3.2) 



of (3.1 1 j in which a%j and &y are constants that influence the relative orientations of molecules. If = bij — 
0, then pairs of molecules prefer to adopt orientations in which their directors are perpendicular to the line 
through their centers. If and/or bij is not zero, then pairs of molecules prefer to adopt orientations that 
are not perpendicular to the line through their centers. For interaction energies of the form (3.2 1, (A. 45 1 
and ( |A.52| yield 



8(B + C) 



[W?(x) Pll (r 2 )(a u _ bii) _ w 2 (x)p 22 (r 2 )(a 22 - b 22 )]r 3 dr, 



(3.3) 



which shows that the spontaneous mean-curvature H Q generally differs from zero when an ^ bn for i = 1 or 
i = 2. 

For another more complicated proposed interaction energy see the work of Yuan (2010). 



4 Possible further developments 

There are numerous ways in which the results of the previous two sections might be generalized. These are 
mainly related to the weakening of the assumptions made in Section [2] Namely, we will consider possible 
ways to weaken Assumptions (ii) and (iii). 

Let us begin with Assumption (iii), which concerns the orientation of the lipid molecules. Although 
Hclfrich (1973) argues that any energetic contribution due to molecular tilt should be negligible and the 
lipid molecules making up a lipid bilayer are almost ubiquitously depicted as being perpendicular to the 
surface formed by the bilayer (see, for example, McMahon & Gallop (2005)), there is no reason to rule out 
the possibility that molecules may tilt relative to the surface normal. Moreover, any such tilt need not be 
small. For simplicity, consider only leaflet 1 and, for brevity, drop the subscript 1. Analogous statements 
apply to leaflet 2. A possible approach to modeling a lipid monolayer made up of tilted molecules is to 
assume that the orientations of those molecules on the surface to be given by a vector field d, so that d(x) 
gives the orientation of the molecules at the point x on the surface representing the monolayer. This vector 
field should satisfy 

|d| = 1 and d • n > 1. (4.1) 



While the first condition in (4.1 1 constrains the lipid molecules to be of fixed length, the second condition 



ensures that the lipid molecules point roughly in the same direction as n. If (4.1 \ 2 did not hold, then either 
the lipid molecules lie flat on the surface (d • n = 0), which seems physically unrealistic, or the tail of the 
lipid molecules are in contact with the adjacent solution (d • n < 0), which, due to the chemical properties of 
the lipid molecules, is physically unattainable. Recall that while the tail of a lipid molecule is hydrophobic, 



G 



the head is hydrophilic. The orientation d should take the place of n in expressions such as (2.6 1. It is 
possible to find a second-order tensor field K such that 



d = Kn. 



(4.2) 



In essence, the tensor K describes the rotation of n into d. Notice that the condition (4.2 1 is not sufficient 
to uniquely determine K, however, if the natural condition 



K(y)z = z for all 



z e T y S, 



yeS 



(4.3) 



is also imposed, where T y S is the tangent space to S at y, then K is uniquely determined. With the 
substitution (4.2 1, an argument analogous to that employed in Section [2] should go through. In particular, 
we conjecture that an expression for the free-energy density with the same form as (2.111 should emerge 
but that the corresponding counterparts of the abbreviations (A. 45 |-( [A.48| will involve K and its surface 
gradient at x. 

To also weaken Assumption (ii), namely the assumption that the lipid molecules are inextensible, two 
steps are necessary. First, the condition (4.1 )i must be removed. Second, the potential governing the 



interaction between the molecules in (4.4) must be replaced by 



$ d (x - y, d(x), e(y)) = ^e~ 2 \x ~ y\\ (x - y) ■ d(x), (x - y) ■ e(x),d(x) ■ e(x), \d(x)\ 2 , \e(y)\ 2 ). (4.4) 

Once again, an analysis similar to that appearing in Section [2] should be feasible and we believe a result 
similar to (2.111 might emerge from such an analysis. 



It should also be feasible to relax Assumption (iv) by adding to the right-hand side of (2.10 1 a term that 
accounts for interactions between the lipid bilayer and the ambient solution. In pursuing such a generaliza- 
tion, it might make sense to follow the lead of Paunov et al. (2000). 

Steigmann (1999) developed a general theory of lipid bilayers, regarded as fluid surfaces with bending 
elasticity. Aside from accommodating a completely general dependence of the free-energy density on the 
mean and Gaussian curvatures, Steigmann's (1999) theory accounts for the areal stretch of the lipid bilayer 
relative to some reference placement. It would be interesting to see if it would be possible to use an argument 
like that contained in Section [2] to obtain a free-energy density that incorporates dependence on the areal 
stretch. To achieve this, it would perhaps be necessary to introduce a reference placement for the lipid 
molecules and then to account for interactions relative to that reference placement. Along with the other 
ideas mentioned in this section, we leave the details of any such analysis to a subsequent work. 



5 Summary 

In Section [2] we presented a microphysical derivation for the Canham-Hclfrich form of the free-energy den- 
sity of a lipid bilayer that includes not only spontaneous mean-curvature but also spontaneous Gaussian- 
curvature. Our approach was inspired by the work of Keller & Merchant (1991). The resulting expressions 
for the spontaneous curvatures allowed us to remark in Section [3] that if the two leaflets are identical and 
the lipid molecules are evenly distributed across the surface of the bilayer, then the natural shape of that 
surface is flat. Moreover, in Section [4] we observed that it should be feasible to extend our method to allow 
for molecular tilt, molecular extension, interactions between the bilayer and the adjacent solution, and areal 
stretch. 



A Appendices 

In these appendices we present the necessary linear algebra and differential geometry necessary to carry out 
the expansions needed in Section [2] We also include some details of the expansion. 

A.l Linear algebra 

Here we recall some basic properties of tensors. Let V be a finite-dimensional vector space equipped with an 
inner-product. For the purposes of this work, a tensor T of order m is a multilinear mapping that takes in 
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m vectors and gives a real number. A vector can be viewed as a first-order tensor and vice versa. A linear 
mapping from V to itself can be viewed as a second-order tensor and vice versa; we denote the set of all such 
linear mappings by Lin(V, V). Other such statements regarding higher-order tensors can also be made. For 
example, a fourth-order tensor may be viewed as a bilinear mapping from V to the space of second-order 
tensors and vice versa. 

Given a tensor T of order m and a tensor S of order n, it is possible to form their tensor product T ® S, 
which is a tensor of order m + n defined by 

T(g> S(vi, . . . ,v m ,v m+ i, . . . ,v m+ „) 

= T(vi,...,v m )S(v m+ i,...,v m+n ) for all v 4 € V, i € {1, . . . , m + n}. (A.l) 

The tensor product is associative, meaning that (T ® S) ® R = T ® (S ® R) and, thus, the use of parenthesis 
in this context may be dropped. 

Just as it is possible to take the inner-product of two vectors, it is possible to take the inner-product of 
two tensors of the same order. The details of how this can be done in general will not be spelled out here. 
However, there are two useful identities that we will need regarding the inner-product of tensors. Namely, 
for any vector z and second-order tensor L, the following hold: 

(z <g> z) ■ L = z • Lz (A. 2) 

(z ® z ® z ® z) • (L ® L) = (z • Lz) 2 (A. 3) 

Assume the dimension of V is 2. The following lemma, which can be established using polar coordinates 
and the integral gradient theorem, involving second- and fourth-order tensors will also be needed. Here, 1 
denotes the identity mapping on V and, thus, is a second-order tensor. 

Lemma A.l. Let g : P — > K be a continuous function. Then 



/ #(|z|)z (g) zda z = ( 7T / g(r)r 3 dr\l, 

J\x\<t \ JO J 



(A.4) 



[ ffflzQz ® z ® z ® zda z = [ g(r)r 5 dr)l ® 1. (A.5) 

J|z|<<? V 4 JO / 

A. 2 Geometry of surfaces 

Consider a surface 5 in a three-dimensional Euclidean point space £ with associated vector space V. Given 
a point j/ in S, let 7^,5 denote the tangent space of iS at y. Let n denote a mapping that determines a 
unit-normal to the surface at each point. Given a mapping h : S — > W defined on the surface that takes 
values in some vector space W, the surface gradient V h of h can be defined by 

Vfft:= V a /i e (l-n(aO®n(x)) for all x G 5, (A.6) 

where /i e is an extension of /i to a neighborhood of x and V^/i 6 is the classical three-dimensional gradient 
of this extension at x. It can be shown that the definition of the surface gradient is independent of the 



extension used on the right-hand side of (A.6 1 



Of particular interest is the opposite L := — V 5 n of the surface gradient V 5 n of n — called the curvature 
tensor, which is a second-order tensor field defined on S. The curvature tensor is symmetric and has two 
scalar invariants: the mean curvature H and Gaussian curvature K defined by 

ff:=itrL, (A.7) 

K :=i[(trL) 2 -tr(L 2 )]. (A.8) 
If K\ and K2 are the two nontrivial eigenvalues of L, often called the principle curvatures, then 

H= 1 -{k 1 + k 2 ) 1 (A.9) 
K = Ki k 2 . (A.10) 



<s 



Figure 2: Depiction of a local curving for a surface S at x. 



A useful identity involving the curvature tensor is provided by the following result. 
Lemma A. 2. For z in V and x in S , 

n(x) ■ (VfLz)z = |L(x)z| 2 . 
Proof. Using the definition and symmetry of L, we have 

n(x) ■ h(x)z = for all x E S. 



(A.ll) 



(A.12) 



Taking the surface gradient of both sides of (A.12 1 in the direction z, using the product rule, and evaluating 
the result at x yields 

- L(.t)z ■ L(x)z + n(x) ■ (Vf Lz)z = 0, (A.13) 



which is (A.ll I. 



□ 



From here on fix a point x in <S. It is possible to parameterize a small neighborhood of S at x using a 
local curving /. Put T := T X S and J\f := (TxS)- 1 , and let U be an open neighborhood of in V. A local 
curving / : hi — > £ of S at x has the following properties: 

1. f(0)=x, 

2. /(u) € S for all u G TnU, 

3. /(u) - + u) G A/" for all u G TflW, 

4. /(u + w) = /(u) + w for all u e T (1W and wejVntf, 

5. V / = 1. 

The gradients of / at describe the local shape of the surface at x. Consider the mappings 



defined by 



F 

A(u,v) 

r(u, v)w 



F : U — ► Lin(V,V), 

A : V x V — > V, 

r : V x V — -> Lin(V,V), 



= v/, 

= [(V F)u]v for all u,v G V, 

= ([(V VF)u]v)w for all u,v,w G V. 



(A. 14) 
(A.15) 
(A.16) 



(A. 17) 
(A.18) 
(A.19) 



The mapping F is a second-order tensor field while the quantities A and T can be viewed as third- and 
fourth-order tensors, respectively. Using Item [3] it can be shown that 



A(zi,z 2 )GA/ r and r(zi,z 2 )z 3 £j\f for all zi,z 2 ,z 3 gT; 



(A.20) 
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further, using Item [4] it can be shown that 

A(zx,n)=0 and r(zi,Z2)n = for all Z],z 2 £ T. 



It follows from ( A.20 1-( A.21 1 that 



trAzi=0 and tr(r(zi, Z2)) = for all Zi,Z2 £ 7". 
The curvature tensor L(x) is related to A by 

A(zi,z 2 ) = n(zi • L(z)z 2 ) for all z 1 ,z 2 eT. 



It follows from |A21| and ( |A.23| that 

Az = n (g) (L(x)z) for all z G T. 

A. 3 Useful approximations 

In this subsection we delineate the key expansions that will be needed in the asymptotic 
energy. The three basic expansions, which follow from the results stated in Appendix |A.2| 

e 2 

/(ez) = x + ez + — (z • L(x)z)n(a;) + o(e 2 ), 
n(/(ez)) = n - eL(x)z - y [(Vf L)z]z + o(e 2 ), 

e 2 

F(ez) = 1 + en(x) ® L(x)z + ~r(z, z) + o(e 2 ). 



Using (A. 11 1 and (A. 25 1-( A.26 1, it transpires that 



e~ 2 \x - f(ez)\ 2 = |z| 2 + -(z • L(x)z) 2 + o(e 2 ), 



(x - /(ez)) • n(x) = - -(z • L(x)z) + o(e 2 ), 
(a; - /(ez)) • n(/(ez)) = ^(z • L(z)z) + o( £ 2 ), 
n(/(ez)).n( 2 ;) = l-^|L(x)z| 2 + o(e 2 ). 



It can be shown that the cofactor F c of F is given by (Gurtin et al. (2010)) 

F c := [F 2 - (trF)F + ^[(trF) 2 - tr(F 2 )]l] T . 



Using (A.lll, (A.22I, ( A.26 1-( A.27 1, and (A.32I, it follows that 



n(/(ez)) • F c (ez)n(x) - 1 + ^|L(x)z| 2 + o(e 2 ). 
Given a scalar-valued function W defined only on the surface, if we write 

W'(x):=V%W and W"(x) := Vf V 5 TU, 

then 



(A.21) 
(A.22) 
(A.23) 
(A.24) 

expansion of the 
are 

(A.25) 
(A.26) 
(A.27) 

(A.28) 
(A.29) 
(A.30) 
(A.31) 

(A.32) 

(A.33) 
(A.34) 



TU(/(ez)) = W(x) + eW'(x)z + — z • W"(x)z + o(e 2 ). 



(A.35) 
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A. 4 Expansion of the energy 



Here we use the results of the previous two subsections to perform the expansion of (2.101 in powers of e. 
Since the domains of integration of the integrals (2.6 1(2.9 1 depend on e through d, this cannot presently be 
achieved. To proceed, we introduce a change of variables that transfers this dependence to the integrand. 



To achieve this, we parameterize Sd{x) using the notion of a local curving / introduced in Appendix A. 2 
which is made possible by (2.2 1. In particular, it is convenient to introduce the set 



75(0) := {z G T X S : |z| < £}, 



(A.36) 



where T X S denotes the tangent space of S at x. Focusing on ( |2.6| and using using the abbreviation ( A.17| ), 
the change of variables y — /(ez) yields 

^u(ar) - e- 2 f <S> d n (x - f(ez),n(x),n(f(ez)))W 1 (x)W 1 (f(ez)) ■ n(/(ez)) • F c (ez)n(x) da z . (A.37) 

JTAO) 



'7H0) 

Let us focus our attention on the quantity 
<(x-/(ez),n(:r),n(/( e z)) = 



cf>n(e- 2 \x - f(ez)\ 2 , (x - /(ez)) • n(x), (x - /(ez)) • n(/(ez)), n(a:) • n(/(ez))). (A.38) 



For all < s < £, put 
and 



^ii,fc(s) := 



0ii,o(«) :=0n(s 2 , 0,0,1), 

01l(6,6,6,6) 



96 



fee {1,2,3,4}. 



(A.39) 
(A.40) 



(6,«2,«3,«4) = (s 2 ,0,0,l) 



Using the expansions ( A. 28 H A. 31 1 and the notation ( A.39 )-( A.40 1, we have 
$Ux - /(ez), n(x), n(/(ez)) - </>i M (|z|) 



+ 



e 2 /l 



(-0 11 ,i(|z|)(z ■ L(x)z) 2 - 0n, 2 (|z|)(z • L(x)z) 

+ 0ii )3 (|z|)(z • L(x)z) + 114 (|z|)|L(a;)z| 2 ) + o(e 2 ). (A.41) 



where L = — V 5 n denotes the curvature tensor of S (see Appendix A.2|. Combining this expansion with 
(A. 33 1 and using the expansion (A. 35) for W\ yields an expansion for the right-hand side of (A.37 1: 



<M*) = f UiiAM)w?(x) + e -Wi(x)[w 1 ( x )(ci> 11 M) 

•'T x 5i(0) V Z 



fa 1A (M))\Hx)z\ 2 + 0n,o(|z|)z • W^'(z)z + - W^foiAMK* ' L W Z ) 2 

+ Wi(a)(0ii,3(|a|) " ^n, 2 (|z|))(z • L(x)z)]) da z + o(e 2 ). (A.42) 
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Similarly, the right-hand sides of (2.7 1— (2.9| can be expanded in powers of e, using notation similar to that 



T e (0) 



of (|A.39|-(|A.40|. Putting these expansions together, we obtain 

[<hiA\'\) w i(x) + 022,o(M)W2 (as) + 2<f )12fi (\z\)W 1 (x)W 2 (x)] da z 

T f (0) 

^22,a(r))](z • L(x)z)da z 

[W?(aO(</>ii,o(|z|) - 0n. 4 (|z|)) + W^oOOWW) 

T f (0) 

022, 4 (|z|)) + 2W 1 {x)W 2 {x){q> 12fi (\z\) + 0i2, 4 (|z|))] |L(x)z| 2 da z 



T f (0) 



[W?(s)0i M (M) + l^(aO0a2,i(M) 



+ 2W 1 (z)W 2 (z)</'i2,i(|z|)] (z ■ L(a;)z) 2 da z 



[(Wi(aO&i j0 (W) + ^ 2 (x)0 12iO (|z|))(z • <(.t)z) 
{W 2 {x)^ 22 .o + Wi(x)</» 12i0 (|z|))(z • W^'(a:)z)] da z . 



(A.43) 



Notice that the first term on the right-hand side of (A.43) is independent of e, while the others are propor 



tional to e 2 . Of the terms involving e, the first is linear is L and the second and third are quadratic in L. 



The last term in (A.43 1 is independent of L. Using the identities (A.2|-(A.5l, the abbreviations 



i> {x) :=2ir / |>ii,o(r)W?(x) + fe.oW^fi) + 2(j> 12 , (r)W 1 (x)W 2 (x)]r dr, 



A(x) 
B(x) 

C{x) 
D(x) 



2 Jo 



[W?(a:)(^u,3(r) - 0ii, 2 (r)) - W 2 (x)(fo 2 , 3 (r) - 22 , 2 (r))]r 3 dr, 
[W?(x)(0 u ,o(r) - <f>u A (r)) + W 2 {x){cp 22fi {r) - 022,4 (r)) 



+ 2W x {x)W 2 (x){cf> l2fi {r) + ^i 2 , 4 (r))]r 3 dr, 
3tt 



1G 



2 Jo 



[Wl 2 (*)0ll,l(r) + ^ 2 2 (x)(/) 22 . 1 (r) + 2W 1 (x)W 2 (x)q> 12A {r)]r 5 dr, 
[(Wi(a;)^ii,o(T-) + W 2 (x)ct> 12t0 (r))tr{W{'(x)) 



+ (W 2 (x) ( t >22 , + W 1 (x)(t>i 2 Ar)MW^(x))]r 3 dr : 
and motivated by the grouping in ( |A.43[ ), we find (on suppressing explicit dependence on x) that 

ip = ip + [AtrL + £tr(L 2 ) + C(trL) 2 + D]e 2 +o(e 2 ). 



Further, upon using (A.7l-(A.8) and the additional abbreviations 

k := 8e 2 (B + C), 
R := ~2e 2 B, 



H := 
K n := 



.4 



A{B + cy 

4D(B + C) - A 2 
8B(B + C) : 



(A.44) 
(A.45) 

(A.46) 
(A.47) 

(A.48) 

(A.49) 

(A.50) 
(A.51) 

(A.52) 
(A.53) 
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allows us to write the expansion of the right-hand side of ( 2.10| a^ 

1 



ip = Tp + ^(H-H o ) 2 + R(K-K o ), (A.54) 
where the correction of o(e 2 ) has been neglected. 
Acknowledgement 

We thank Mohsen Maleki for very fruitful discussions. 

References 

[1] Canham PB (1970) The minimum energy of bending as a possible explanation of the biconcave shape 
of the human red blood cell. Journal of Theoretical Biology 26:61-81 

Cosserat E, Cosserat F (1909) Theorie des Corps Deformables. Herman et fils, Paris 

Dobereiner HG, Selchow O, Lipowsky R (1999) Spontaneous curvature of fluid vesicles induced by 
trans-bilayer sugar asymmetry. European Biophysics Journal 28: 174-178 

Foppl A (1907) Vorlesungen iiber technische Mechanik, Bd. 5, Die wichtigsten Lehren der hoheren 
Elastizitatstheorie, Teubner, Leipzig 

Germain S (1821) Recherches sur la Theorie des Surfaces Elastique. Huzard-Courcier, Paris 

Gurtin ME, Fried E, Anand L (2010) The Mechanics and Thermodynamics of Continua. Cambridge 
University Press, New York 

Hclfrich W (1973) Elastic properties of lipid bilayers: Theory and possible experiments. Zeitschrift fur 
Naturforschung 28c:693-703 

Keller JB, Merchant GJ (1991) Flexural rigidity of a liquid surface. Journal of Statistical Physics 
63:1039-1051 

Lasic DD (1988) The mechanism of liposome formation. A review. Biochemical Journal 256:1-11 

Lipowsky R (1990) Shape fluctuations and critical phenomena. In van Beijeren, H (ed) Fundamental 
Problems in Statistical Mechanics VII. North-Holland, Amsterdam, pp 139-170 

Ljunggren S, Eriksson JC (1985) Comments on the origin of the curvature elasticity of vesicle bilayers. 
Journal of Colloid and Interface Science 107:138-145 

Luisi PL, Walade P (2000) Giant Vesicles. Wiley, Chichester 

Mai'er W, Saupe A (1958) Eine einfache molekulare Theorie des nematischen Kristallinfliissigen Zus- 
tands. Zeitschrift fur Naturforschg 13a:564-566 

Maleki M, Seguin B, Fried E (2012) Kinematics, material symmetry, and energy densities for lipid 
bilayers with spontaneous curvature. Biomechanics and Modelling in Mechanobiology 10.1007/sl0237- 
012-0459-7 

McMahon HT, Gallop JL (2005) Membrane curvature and mechanisms of dynamic cell membrane 
remodelling. Nature 438:590-596 

Paunov VN, Sandler SI, Kaler EW (2000) A simple molecular model for the spontaneous curvature and 
the bending constants of nonionic surfactant monolayers at the oil/water interface. Langmuir 16:8917- 
8925 



1 See d A. 7b and 1A.81 for the definitions of the mean H and Gaussian K curvatures. 



13 



[17] Poisson SD (1812) Memoire sur les surfaces elastiques. Memoire de Classe des Sciences Mathematiques 
et Physiques de I'Institut de France 2nd pt.: 167-225 

[18] Sackmman E (1994) The seventh Datta Lecture. Membrane bending energy concept of vesicle- and 
cell-shapes and shape-transitions. FEBS Letters 346:3-16 

[19] Safran SA (1994) Statistical Thermodynamics of Surfaces, Interfaces, and Membranes. Addison-Wesley, 
Reading 

[20] Seguin B, Fried E (2012) Statistical foundations of liquid-crystal theory I. Discrete systems of rod-like 
molecules. Archive for Rational Mechanics and Analysis 206:1039-1072 

[21] Seifert U (1997) Configurations of fluid membranes and vesicles. Advances in Physics 46:13-137 

[22] Steigmann DJ (1999) Fluid films with curvature elasticity. Archive for Rational Mechanics and Analysis 
150:127-152 

[23] von Kaman, Th (1910) Festigkeitsprobleme im Maschinenbau. In: Klein F, Miiller C (eds) Encyklopadia 
der mathematischen Wissenschaften IV/4, Teubner, Berlin, pp 311-385 

[24] Winterhalter M, Helfrich W (1988) Effect of surface charge on the curvature elasticity of membranes. 
Journal of Physical Chemistry 92:6865-6867 

[25] Yuan H (2010) A solvent-free coarse-grained model for biological and biomimetic fluid membranes. Ph. 
D. Thesis. Pennsylvania State University 



14 



